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ABSTRACT 

We report the results of a study exploring the stellar populations of 13 luminous (L 
CN > 1.2L*), spectroscopically confirmed, galaxies in the redshift interval 5.5 < z < 6.5, 

all with Hubble Space Telescope (HST) WFC3/1R and SpitzerlRAC imaging from the 
HST/CANDELS and Spitzer/SEDS surveys. Based on fitting the observed photometry 

j» with galaxy spectral energy distribution (SED) templates covering a wide range of 

^si different star-formation histories, including exponentially increasing star-formation 

rates and a self consistent treatment of Lyman-alpha emission, we find that the derived 
stellar masses lie within the range 10 9 M Q < M* < 10 10 M Q and are robust to within 
a factor of two. In contrast, we confirm previous reports that the ages of the stellar 
populations are poorly constrained. Although the best-fitting models for 3/13 of the 
. . sample have ages of > 300 Myr, the degeneracies introduced by dust extinction mean 

that only two of these objects actually require a > 300 Myr old stellar population 
to reproduce the observed photometry. Moreover, when considering only smoothly- 
varying star- format ion histories, we observe a clear tension between the data and 
models such that a galaxy SED template with an old age is often chosen in order 
to try and fit objects with blue UV-slopes but red UV-to-optical colours. To break 
this tension we explore SED fitting with two-component models (burst plus on-going 
star- formation) , thereby relaxing the requirement that the current star formation rate 
and assembled stellar mass must be coupled, and allow for nebular line+continuum 
emission. On average, the inclusion of nebular emission leads to lower stellar-mass 
estimates (median offset 0.18 dex), moderately higher specific star- format ion rates, 
and allows for a wider range of plausible stellar ages. However, based on our SED 
modelling, we find no strong evidence for extremely young ages in our sample (i.e. 
< 50 Myr). Finally, considering all of the different star-formation histories explored, 
we find that the median best-fitting ages are of the order ~ 200 — 300 Myr and that 
the objects with the tightest constraints indicate ages in the range 50-200 Myr. 

Key words: galaxies: high-redshift - galaxies: evolution - galaxies: formation. 
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1 INTRODUCTION 



A combination of near-Infrared and Spttzer (Werner et al 



2004| ) Infrared Array Camera (IRAC; |Fazio et al.| |2004 1 



observations have shown that many high-redshift (z ^ 
5) Lyman-break galaxies (LBGs) have red rest-frame 
ultraviolet-optical colours (e.g. Yan et al.|2005 Eyles et al 



2005| |Ono et al.||2010| |Labbe et al.|2010[ ). If this red colour 
is taken to indicate the strength of the Balmer break, then 
old stellar ages (> 300 Myr) are inferred for these objects, 
indicating a formation redshift at z > 8. This property was 
first reported for two i— drop galaxies in |Eyles et al.| |2005), 
with their spectral energy distribution (SED) fitting indi- 
cating ages in the range 250-650 Myr, and stellar masses of 
order 20% of current-day L* galaxies. A later analysis of a 
larger sample showed that 40% of their sample display ed ev- 
idence for substantial 4000A/Balmer spectral breaks (Eyles 
et al.||2007| ). More recently, |Richard etH] \2Qll\ reported 
the discovery of a lensed galaxy at z = 6.027 which also 
appears to have a strong Balmer-break, consistent with a 
mature stellar population. If confirmed, a substantial popu- 
lation of galaxies at z > 6 with mature stellar populations 
has important consequences for star formation at z > 8 and 
its contribution to reionization. 

SED-fitting techniques, however, often employ tem- 



plates that only account for stellar emission (eg. |Bruzual fc| 
Chariot 2003 Maraston 2005 I and recently several groups 



have investigated how the inclusion of nebular emission in 
the fitting can affect the derived parameters (e.g. IRobert- 



son et al.||20T0l |Schaerer fc de Barros||2010| |Labbe et al 



2010||Ono et al.|2010[ ). For example, |Labbe et al.| ( |2010[ ) re^ 
ported difficulties in reconciling the observed Balmer break 
with a very blue ultraviolet (UV) continuum slope (f3 ~ —3; 
fx oc A* 9 ) observed in a stack of the photometry for the 
faintest galaxies in a sample of z ~ 7 LBGs (Hiso.ab > 27.5) 
taken from the Hubble Ultra-deep Field (HUDF) and Early 
Release Science (ERS) field. They suggest that episodic star 
formation may be required to match the colours of these ob- 
jects (being simultaneously blue in the rest-frame UV and 
red in colours spanning the 4000A/Balmer break), as well 
as the contribution of nebular emission lines that are not 
included in the standard population synthesis codes. An 
illustration of the potentially dramatic effect that nebular 
emission can have on derived stellar masses and ages is pro- 



vided by Ono et al. (2010), in their fitting of stacked multi- 



wavelength observations of Lyman alpha emitters (LAEs) 
at z ~ 5.7 and z ~ 6.6 within the Subaru/ XMM-Newton 
Deep Survey (SXDS) field. The derived masses from mod- 
els with maximum nebular emission (escape fraction of 
ionising photons, f esc = 0), are more than an order-of- 
magnitude lower than those based on only stellar emission 
(e.g. M„ ~ 3 x 10 7 M compared to M* ~ 5 x I0 8 M for the 
z ~ 5.7 stack) and the derived ages are also much younger 
(~ 3 Myr compared to ~ 300 Myr). 

Recent results from cosmological hydrodynamic simu- 
lations of reionisation predict that Lyman-break selected 
galaxies have young ages, low intrinsic reddening and sub- 
solar metallicities at high redshifts (e.g. ~ 50 — 150 Myr 
at 6 < z < 8; Finlator et al. 2010 I with the distribution 



broadening to older ages at later times (e.g. 200-600 Myr at 
Finlator et al.||2006[ ) . The results from |Dayal et al.| 



logical smoothed-particle hydrodynamic (SPH) simulation, 
also support these predictions. By matching models to ob- 
served Lya and UV-luminosity functions at 5.7 < z < 7.6, 
the brightest LAEs, which are shown in their simulation to 
be a subset of the LBG population ( Dayal fc Ferrara|2011 1, 
are expected to have intermediate ages (~ 200 Myr), with 
the spread in age increasing to lower luminosities (ages from 
a few Myr to 300 Myr). 

Given that it is difficult to constrain the star-formation 
history (SFH) of high-redshift galaxies from the available 
photometric data, many studies investigating the physical 
properties of high-redshift galaxies use small template sets 
(including restricted ranges in reddening and metallicity) 
when SED fitting to avoid introducing too many degenera- 
cies. In this way, if the template sets are consistent, the 
results for samples at different redshifts can be compared 
to one another to search for any evolution in the typical 



galaxy properties. The study by Yabe et al. (2009) illus- 



trates the importance of matching the template sets when 
comparing the physical properties of their z ~ 5 LBGs to 
similar samples at different redshifts, where differences be- 
tween the template sets affect the distributions of derived 
star-formation rates (SFRs) and ages by up to a factor of 
10. They do show, however, that derived masses are fairly 
insensitive to changes in template SFHs. 

This approach is limited in its usefulness, however, if the 
galaxies evolve in a parameter space that is not allowed to 
vary in the template set, such as in dust content, metallicity 
or even in the SFHs. Another limitation is in the search for 
trends in physical properties within the individual samples 
(eg. SFR vs. stellar mass), as any apparent trends may be 
imposed by the template set itself, as discussed by |McLure| 



et al. (2011 1. When comparing to other galaxy populations, 



for example to determine how high-redshift LBGs are linked 
to lower redshift populations, it is clearly desirable to have 
robust parameter determination that is independent of the 
adopted fitting method. 

Tracking the evolution of these trends gives us the 
best indication of the typical mass-assembly route for these 
galaxies. In particular, tracking the SFR (or UV luminos- 
ity) of galaxies of a given mass over time can rule out cer- 
tain SFHs. Stark et al. 12009 ) tracked the evolution of the 
physical properties of LBGs from z ~ 6 to z ~ 4, find- 
ing little evolution in the stellar mass-UV luminosity rela- 
tion (M* — Luv) over this redshift range. They suggested 
that this is best explained by the different redshift bins be- 
ing populated by distinct populations with short periods of 
star formation, as opposed to the population remaining in 
the LBG phase, constantly forming stars and building up in 
stellar mass over the entire redshift range. This view is also 
supported by the clustering analysis of star-forming galaxies 
performed by Lee et al. (20091, which suggested a typical 
duration of star formation < 400 Myr (la) at z ~ 4 and 
z ~ 5. 

The alternative explanation, however, is that these 
galaxies exhibit rising SFHs, a feature predicted by the mod- 



els of Finlator et al. (20101 and supported by the obser- 



1 2009 ) , in which Lya emission is incorporated into a cosmo- 



vation that high-redshift galaxies selected to have constant 
number density display UV luminosities which increase with 
time ( jPapovich et al.||2010 ). Moreover, an exponentially in- 
creasing SFR would naturally produce little evolution in the 
observed specific SFR (sSFR), which current data suggests 
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remains fairly constant at ~2.5 Gyr -1 from z ~ 2.5 — 7 
( Gonzalez et al.||2011 McLure et al.pbll l. However, these 
observations provide a tension with current models of galaxy 
formation and evolution which predict sSFRs that continue 
to rise with increasing redshift proportional to (1 + z) 2 ' 25 



2011) 



(Neistein & Dekel 2008 Dave, Oppenheimer & Finlator 



A different approach to determining the physical prop- 



erties by SED fitting is used by Schaerer & de Barros (2010), 



who allow for a large range of SFHs, dust extinction, metal- 
licity and nebular emission. They show that a wide range of 
ages and SFRs are acceptable within the confidence contours 
of their fits, although they also have the added complications 
of uncertainty in photometric redshift and Lya flux inherent 
to photometrically-selected dropout samples compiled from 
the literature. 

The new Hubble Space Telescope (HST) WFC3 data 
taken as part of the Cosmic Assembly Near-infrared (IR) 



Deep Surve y (CANDELS, |Grogin et aL]|2011| |Koekemoer 
|et al.||2~011[ ), combined with mid-IR data from the Spitzer 
Extended Deep Survey (SEDS) allows us to place tighter 
constraints on the rest-frame UV and optical fluxes of high- 
redshift galaxies. Motivated by the newly-available data, in 
this paper we explore the potential to accurately constrain 
the physical properties of high-redshift galaxies based on a 
sample of spectroscopically confirmed LBGs in the redshift 
range 5.5 < z < 6.5. 

The structure of the paper is as follows. In section 2 we 
describe the sample as well as the determination of the Lya 
equivalent widths (EWs) from the available spectroscopy. In 
Section 3 we outline the photometric data used to construct 
the sample SEDs, as well as the deconfusion of the IRAC im- 
ages. In Section 4 the observed characteristics of the objects 
are examined and compared to a more complete photomet- 
rically selected sample. The SED fitting method is outlined 
in Section 5 along with a summary of the different tem- 
plate sets used. Section 6 gives the results from fitting with 
the smoothly varying SFHs starting from more restricted 
template sets to compare results directly to the literature, 
then with a wider range of models and parameters. Section 

7 details the results of fitting with simple two-component 
templates both with and without nebular emission. Section 

8 provides a discussion of the derived physical parameters 
and our conclusions are presented in Section 9. Through- 
out this paper we assume a cosmology with Ho = 70 km 
s^Mpc^ 1 , Q m = 0.3, fl A = 0.7. All magnitudes are quoted 
in the AB system ( |Oke fc Gunn||1983 1. 



2 SPECTROSCOPIC DATA 

2.1 Spectroscopic samples in the UDS and 
GOODS-S fields 

The properties of the 13 LBGs used in this paper are sum- 
marised in TableQ] The fundamental criteria for inclusion in 
the sample were that LBG-selected objects had z spec > 5.5 
and m z < 26.5 and were covered by both WFC3 near-IR 
and deep Spitzer 3.6//m and 4.5/xm imaging. Objects were 
drawn from separate two sources. 

Eleven of the objects are taken from the spectroscopic 



one of the five fields covered by CANDELS. This particular 
study is chosen because it currently provides the most com- 
plete sample of spectroscopically targeted LBGs at z ~ 6 
with spectra of high enough quality to determine Lya prop- 
erties. 



The objects were selected by requiring 



> 5.5, 



m z < 26.5 as well as requiring > 5<r detections in the WFC3 
H— band and a detection in both the IRAC 3.6/im and 
4.5pm bands. The first two criteria allowed us to select UV- 
luminous LBGs in a narrow redshift range while the IRAC 
detections ensured that the SED could be sampled red- 
wards of the 4000 A/Balmer- break, essential for constraining 
masses. This selection left us with 4 objects selected from 
the Veo6-band dropout sample, and 7 from the i775-band 
dropouts. One object meeting the above requirements was 
rejected due to it's proximity to a nearby star that would 
make deconfusion of the IRAC bands incredibly uncertain. 
Furthermore, three objects that were bright enough in zgso 
were not bright enough in the H— band for use as reliable 
priors in the IRAC deconfusion, and so were also left out of 
the sample along with one final object that lacked a detec- 
tion at 4.5/im. 

We also required the spectroscopic redshift quality flag 
to be A or B, although after visual inspection, two absorp- 
tion systems were included which had a quality flag of C in 



Vanzella et al. (2009). All but one of the absorption systems 



with z avec > 5.5 were assigned a quality flag of C by |Vanzella| 



et al. ( 2009 1, given the difficulty of unambiguously assigning 



the observed continuum break to Lya absorption. However, 
in the interest of examining the properties of a sample of 
z ~ 6 objects covering a full span in intrinsic properties, we 
have kept them in our sample. We note that the high quality 
WFC3 near-infrared plus Spitzer IRAC data provide pho- 
tometric redshifts for these objects which are in excellent 
agreement with the spectroscopic redshifts, indicating that 
they are robust. 

This sample is supplemented by two objects, taken from 
a sample of spectroscopically confirmed z > 6 LBGs in the 
Ultra Deep Survey (UDS) as presented in|Curtis-Lake et al.| 



(20121, which also fall within the area targeted by CAN- 



DELS. The objects' parent sample of LBGs was selected us- 
ing a photometric-redshift analysis (see |McLure et aL||2009| 
for details) and the spectroscopic sample is described in 
Curtis-Lake et all (12012 1. 



2.2 Lya properties 

We measure the Lya EW for each object from the spectra. 
The main aim is to determine a robust estimate that can be 
used to add the Lya flux to the models before SED fitting. 
As such we chose a method that links the EW to the multi- 
wavelength photometry, using an approach similar to that 



described in Curtis-Lake et al. (2012 I, to which the reader is 



sample of Vanzella et al. (2009) within the GOODS-S field, 



referred for a fuller explanation. Essentially, the spectrum is 
first flux-calibrated by integrating over the filter profile that 
encompasses Lya and scaling it to match the photometry 
in that band. To prevent sky-line residuals dominating the 
integrated filter flux, the continuum is modelled by a sim- 
ple step function with the continuum estimated from clean 
regions of the spectrum red-wards of Lya. The error in the 
continuum estimate is folded into the final error on the de- 
rived EW. 
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Table 1. Full sample of objects for SED fitting in the GOODS-S and UDS fields taken from Vanzella et al. (2009) and Curtis-Lake et 
al. (2012) respectively. The first column provides the ID used to identify the objects throughout this work. The second column gives the 
GOODS/UDS ID incorporating the J2000 coordinates. Columns 3-5 give the spectroscopic redshift, object class (emission line (em) or 
absorption (abs) system) and redshift quality flag (A-C, see text for details). Columns 6 and 7 give the total z— band apparent magnitude 
along with the luminosity in units of L* (taken to be (6.3 ± 0.7) X 10 40 erg s — 1 A -1 from the UV luminosity function at z ~ 6, McLure 
et al. 2009). Column 8 gives the paper reference; Vanzella et al. 2009 (V09) or Curtis-Lake et al. 2012 (CL12). The final column gives the 
broad-band calibrated EW (see Section 2.3 for details) for each of the emission line objects where there was visible continuum red-wards 
of Lya in the spectrum. 



Our ID 


GOODS/UDS ID 


Zspec 


Class 


Quality 


m z 


L (/ L*) 


Reference 


EW rest (/A) 


2 


J033215.90-274123.9 


5.57 


cm 


B 


25.84 


1.9 


V09 


6.9 ± 1.1 


5 


J033225.61-275548.7 


5.79 


em 


A 


24.58 


6.4 


V09 


6.5 ± 1.1 


6 


J033246.04-274929.7 


5.79 


em 


A 


25.98 


1.8 


V09 


31.0 ± 4.0 


8 


.1033240:01-274815.0 


5.83 


cm 


A 


25.23 


3.6 


V09 


14.8 ± 2.0 


13 


J033228.19-274818.7 


5.94 


cm 


B 


26.43 


1.2 


V09 


24.8 ± 4.9 


18 


J033224.80-274758.8 


5.99 


cm 


B 


26.10 


1.7 


V09 




23 


J033222.28-275257.2 


6.20 


abs 


C 


26.17 


1.6 


V09 


n/a 


24 


J033224.40-275009.9 


5.50 


abs 


C 


25.67 


2.2 


V09 


n/a 


25 


J033237.63-275022.4 


5.52 


cm 


A 


26.03 


1.6 


V09 


8.2 ± 1.5 


26 


J033211. 93-274157.1 


5.58 


cm 


B 


26.35 


1.2 


V09 




27 


J033245.23-274909.9 


5.58 


cm 


B 


25.94 


1.7 


V09 


31.6 ± 7.0 


157 


UUDS_J021800.90-051137.8 


6.03 


cm 


A 


24.85 


5.3 


CL12 


54.7 ± 11.3 


248 


UUDSJ021735. 34-051032.6 


6.12 


cm 


A 


25.13 


4.2 


CL12 


46.8 ± 8.4 



The continuum at Lya is then estimated using either 



Subaru NB921 narrow-band imaging (Sobral et al. 20111, 



providing photometry just red- wards of Lya for the two 
objects in the UDS, or HST WFC3 F125W imaging for 
the objects in the GOODS-S, assuming a fiat slope in F v 
[F\ oc A~ 2 ) in each case. Although, in principle, Lya EWs 
can be derived directly from spectra without the need to 
flux calibrate, the S/N in the measured continuum is low 
for these faint objects, and linking the measured EW with 
the photometry allows us to reliably add Lya to the models 
in a self-consistent fashion. 



Two of the objects from the Vanzella et al. (2009) sam- 



ple (IDs 18 and 26) show no visible continuum red-ward of 
Lya and as a consequence any attempt at flux calibration is 
very uncertain. For these objects the filters that have Lya 
contributing to the filter flux are not used in the SED fitting. 



3 PHOTOMETRIC DATA 

In this section we provide a brief description of the various 
optical, near-IR and mid-IR imaging datasets which were 
used to construct the SEDs of our sample of z ~ 6 LBGs. 
In all cases the optical and near-IR HST photometry was 
measured using 0.6"— diameter circular apertures, with the 
WFC3/IR data in the F125W(J 125 ) and F160W(H WO ) fil- 
ters corrected for missing flux relative to the optical ACS 
bands. The ground-based photometry was measured in cir- 
cular apertures with diameters of either 2.0" (Subaru and 
UKIDSS) or 1.2" (HAWK-I) depending on seeing (FWHM 
~ 0.8"and ~ 0.5"respectively), in order to match the en- 
closed flux of the Hiao photometry. Obtaining accurate pho- 
tometry from the confused Spitzer IRAC imaging was a 
more complicated process and is discussed further below. A 
full list of the photometry adopted during the SED fitting 
analysis of each individual object can be found in Tables 1 
and 2 of the Appendix. 



3.1 Optical imaging 

For the eleven objects which originate from the Vanzella 
et al. (2009) sample, the optical photometry was mea- 
sured from the publicly available HST ACS imaging of the 
GOODS-S field in the F435W, F606W, F775W and FS50LP 
filters (GOODSv2.0; Giavalisco et al. 2004). For the two ob- 
jects in the UKIDSS UDS field the optical photometry was 
measured from the deep Subaru imaging of the UDS in the 
BVRi'z' filters (Furusawa et al. 2008). 



3.2 Near-IR imaging 

The primary near-IR photometry was measured from the 
publicly-available CANDELS imaging of the GOODS-S and 
UDS fields with WFC3/IR in the J 125 and H 160 filters (Gro- 
gin et al. 2011; Koekemoer et al. 2011). In addition, ground- 
based K— band photometry was measured for the two UDS 
objects using data release 8 (DR8) of the UKIDSS UDS 
and K— band photometry for 7/11 of the GOODS-S objects 
was measured from the first epoch of the HAWK-I imaging 
programme of the CANDELS GOODS-S and UDS fields 
(HANDELS; P.I. A. Fontana). 



3.3 Spitzer IRAC imaging 

A key element of this study is the crucial information on 
the LBG stellar masses supplied by the available Spitzer 
IRAC data. For each of the LBGs the photometry at 3.6pm 
and 4.5pm was measured from either the publicly available 
reductions (v3.0) of the IRAC imaging obtained as part of 
the GOODS survey (proposal ID 194, Dickinson et al., in 
preparation), or the Spitzer Extended Deep Survey (SEDS; 
Ashby et al. 2012, in prep.) imaging of the UDS. 
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3.3.1 Deconfusion analysis 



5 SED FITTING - METHOD 



Due to the depth of the available IRAC imaging, and 
the comparatively broad IRAC PSF (FWHM ~ 1.7"), the 
3.6/mi and 4.5/rni data in the GOODS-S and UDS CAN- 
DELS fields is heavily confused. As a result, it is necessary to 
employ some form of deconfusion analysis to allow accurate 
photometry to be extracted beyond the natural confusion 
limit. 

In this study we employ the deconfusion software devel- 
oped by McLure et al. (2011), which exploits the high spatial 
resolution WFC3/IR -H160 imaging to provide normalized 
templates of each object in the field and then, via a trans- 
fer function, produces a synthetic IRAC image. Through a 
matrix inversion procedure, the code then simultaneously 
fits the amplitude (or total flux) of each template in or- 
der to produce the optimal reproduction of the observed 
IRAC image. Based on the best-fitting model it is then pos- 
sible to extract photometry at 3.6/im and 4.5/im in circular 
0.6"— diameter apertures, matched to the WFC3/IR pho- 
tometry in the H\eo filter. Full details of the deconfusion al- 
gorithm are provided in Appendix A of McLure et al. (2011). 



4 OBSERVED COLOURS AND SAMPLE BIAS 

The available photometry samples three key regions of the 
object SEDs; the rest-frame UV, the rest-frame optical and 
the 4000A/Balmer-break region (which is bracketed by the 
near-IR and Spitzer observations). We plot the information 
we have on these regions in Fig. [I] using observed colours, 
with the observed rest-frame near-UV slope being described 
by J — H, the Balmer-break region (the rest-frame UV-to- 
optical colours) bracketed by H — 3.6/im, and the slope in 
the rest-frame optical by 3.6/im — 4.5/im. The J — H colour 
does not give an accurate measurement of /3 in this context, 
as the filters sample the near-UV red-wards of the region 
from which (3 is estimated (i.e. A rcst ~ 1500A). 

Fig.[l](b) shows that the IRAC 3.6/im — 4.5/im colours 
extend to quite blue colours and that this does not seem to 
be dependent on the colours spanning the Balmer-break. In 
fig. [l](c), however, we see a correlation between J — H and 
3.6/im — 4.5/im. Discussion of the physical interpretation of 
these trends is deferred to Section 8.1. 

Fig. [T] (a) suggests a weak anti-correlation between 
J — H and H — 3.6/im. Also plotted is a more complete 



6 sample (plotted in green, taken from McLure et al 



2011). This LBG sample was selected by a photometric red- 



shift technique within the HUDF and ERS fields and the 
plotted objects are those with 6.0 < z p hot < 6.5. We can see 
that the colours of the underlying sample occupy a similar 
colour-colour space to the majority of our sample, although 
the reddest J — H colours are mostly found in objects which 
also have very red H — 3.6/im colours. A positive correla- 
tion between these colours would be expected if dust were 
reddening the SEDs, and the lack of red H — 3.6/im colours 
in our sample would suggest that the spectroscopically con- 
firmed objects, that in the most part show Lya in emission, 
are relatively dust free. 



We fit a range of star formation histories (SFHs) to the indi- 
vidual object SEDs using the LePhare photometric redshift 
code ( Ilbert et al.|2006 1 . These SFHs include exponentially- 
decreasing star formation rates (r models), constant star for- 
mation, instantaneous burst, exponentially-increasing (EI) 
star formation as well as two component models with an 
old burst plus ongoing constant star formation with and 
without full nebular emission. The different template sets 
employed throughout this paper are su mmarised in Table [2 
and all templates are built using the Bru zual fc Chariot 
(20031 stellar evolutionary models. 



The full set of smoothly varying SFHs (template set C) 
consists of r models (r = 50, 100, 200, 500 Myr, 1, 2, 5, 10 
Gyr), instantaneous burst, constant SFR and exponentially 
increasing (EI) SFR models (with sSFR = 2, 4, 6, 8, 10, 
12, 14, 16, 18, 20 Gyr" 1 at all ages, assuming that the stel- 
lar mass at any time is the integral under the past SFH). 
Ages were allowed to vary between 10 Myr to the age of 
the Universe at the redshift of each object, and metallici- 
ties ranging from solar (Z©) to l/50th solar (0.02 Zq) were 
included. Extinction by dust was applied using the |Calzetti| 
|et al. I ( |2000[ ) dust law for starburst galaxies in 150 steps of 
E(B - V) varying from to 0.5. 

The EI models are made employing an initial seed mass 
with the specific star-formation rate then being constant 
with time. The parameterisation, however, assumes the stel- 
lar mass at any time is the integral under the past star- 
formation history. As the stellar mass, in fact, decreases 
fairly rapidly after a burst of star formation, this assump- 
tion leads to instantaneous sSFRs somewhat higher than the 
parameterised value. Moreover, for low sSFR models with 
young ages, the seed mass can account for a significant frac- 
tion of the mass in the model and the resultant SED is in- 
distinguishable from a large burst followed by constant star 
formation. 

The two restricted template sets (A and B) use subsets 
of the full range of smoothly- varying SFHs to test the effect 
of various assumptions on the derived physical parameters 
and to check consistency with previous studies. In partic- 
ular template set B includes models consistent with those 



adopted by Eyles et al. (2005 20071, who provided some of 



the strongest evidence to date for old stellar populations at 
high redshift. 

The burst+constant star formation rate models (set D) 
were produced in a similar manner to the method applied by 



Eyles et al. ( 2005 2007 1, with the same range of metallicities 



as the smoothly varying SFHs, and star-forming fractions 
varying from 0-15% by mass in stars. This upper limit in the 
mass of the star-forming component was chosen to provide 
a minimum flux contribution to the 3.6/im band of 10%. 
The star-forming fractions were characterised by constant 
star formation of durations 10, 30 and 100 Myr. The age 
of the burst varied between the age of the duration of star- 
formation for that model, and the age of the Universe at 
the redshift of each object. Extinction was applied using a 
coarser grid, with 13 steps in E(B — V) varying from to 
0.495. 

The final template set (E) adds nebular emission, based 
on the Robertson et al. (2010) model, to the burst+constant 



SFR models (D) according to the prescription outlined in 
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Figure 1. The colour-colour space spanned by the sample, in which the rest-frame near-UV slope is described by the J — H colour, 
H — 3.6/im colours span the Balmer-break region and 3.6 — 4.5/im colours describe the slope in the rest-frame optical. The filled diamonds 
are the objects that were originally confused in the IRAC images and the open symbols show the objects that are relatively isolated. A 
photometrically-selected 6.0 < z v ^ ot < 6.5 sample of LBGs taken from McLure et al. (2011) is plotted in green in the left-hand panel. 
Objects from this sample without detections at 3.6/im have their H — 3.6/xm colours plotted as upper limits. 



Table 2. Summary of the different template sets used in the SED fitting. The first column gives the label used in the text to refer to the 
given template set, the second column gives a short description of the template set. Column 3 details the models included and columns 
4-5 give the range of metallicities and extinction explored in each set. The final column displays whether or not the template set includes 
nebular emission. 



Label 


Description 


Templates used 


Metallicities 

(/Z ) 


E(B-V) 


Nebular 
emission (Y/N) 


A 


Restricted template set 


Constant SFR 
Burst 

r = 100, 500 Myr 


0.2 





N 


B 


For comparison to 
Eyles et al. (2007) 


Constant SFR 
Burst 

r = 50, 100, 200, 500 Myr, 1 Gyr 


1 





N 


C 


All smoothly varying SFHs 
plus a burst 


Constant SFR 
Burst 

r = 50, 100, 200, 500 Myr, 1, 2, 5, 10 Gyr 

Exponentially increasing SFR, 

SSFR = 2, 4, 6, 8, 10, 12, 14, 16, IS, 20 Gyr" 1 


0.02, 0.2, 0.4, 1 


150 values from 

- 0.5 
in steps of 0.003 


N 


D 


Burst + recent constant 
star formation 


Underlying burst, mass fraction 

- 15% (0-2% in steps of 0.1%, 3-15% in steps of 1%) 
PLUS recent constant star formation 
with duration = 10, 30, 100 Myr 


0.02, 0.2, 0.4, 1 


13 values from 
- 0.495 
in steps of 0.038 


N 


E 


Burst + recent constant 
star formation with 
nebular emission added 


Underlying burst, mass fraction 

- 15% (0-2% in steps of 0.1, 3-15% in steps of 1) 

PLUS recent constant star formation 


0.02, 0.2, 0.4, 1 


13 values from 

- 0.495 
in steps of 0.038 


Y 



with duration = 10, 30, 100 Myr 



McLure et al. (20111. Intergalactic medium (IGM) absorp- 
tion short-ward of Lya is applied according to the |Madau| 
( 1995 1 prescription, and a Chabrier ( 2003 1 IMF is assumed 



in all cases. 



Lyman-a emission is added to each of the models using 
the measured EW values, using the same filter to provide 
the UV continuum estimate from the models as used for the 
original EW measurements. The Lya emission is added to 
the models after extinction is applied and after IGM absorp- 
tion, to ensure that it is comparable to our measured value. 
It is imperative to add the Lya flux if fitting to the filter 
containing Lya emission because, with the typical EW for 
this sample (~ 25A), it can alter the photometry by 0.2-0.5 
magnitudes depending on the exact redshift. For the two ob- 
jects without reliable Lya EW measurements (18 and 26), 
any filters with a Lya flux contribution are discounted from 
the SED fitting. 



6 SED FITTING - SMOOTHLY VARYING 
STAR FORMATION HISTORIES 
(TEMPLATE SETS A-C) 

6.1 Restricted model set (template set A) 

We begin the process of investigating our ability to constrain 
the physical properties of the high-redshift sample by start- 
ing with the most restricted set of models, as is commonly 
done in the literature (eg. Stark et al.|2009" Gonzalez et al. 
2010 1. This model set is chosen taking the view that there is 



little evidence for much intrinsic reddening in LBGs, as the 
average UV continuum slopes that are observed are quite 
blue ( |Bouwens et al.|2009||McLure et al.|20lT||Dunlop et al 



20121. Consequently, we begin our model fitting by restrict- 
ing our model set to a single sub-solar metallicity (l/5th 
solar) with zero reddening, and consider only constant SFR, 
burst and r models with r = 100 and 500 Myr. 

Table [3] shows the results from fitting with this model 
set, with the errors determined from A% 2 =1, giving 68% 
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Table 3. Results from SED-fitting using a restricted template set (template set A) which features a single metallicity, no dust reddening 
and either constant, burst or r = 100, 500 Myr models. The first column gives the object ID, columns 2-4 give the best-fitting stellar 
mass, SFR and age respectively. The mass and SFRs have been scaled to total using the appropriate H-band aperture corrections. The 
errors quoted are the 68% confidence limits marginalised over all other parameters. The final three columns give the best-fitting model, 
X 2 and the number of filters used in the fit. The typical number of degrees of freedom in the fitting for 9 photometric points is 6. 



ID 




SFR 

(/M yr- 1 ) 


loETinfAtrp /Ovr^ 




Model 


y2 


filters 


2 


a - oo -0.10 


n on+ ' 00 


7.861°;°= 




burst 


3.60 


8 


5 


9 74+0.22 

a - ' -0.08 


14 9Q+ 4 - 85 




r 


= 100 Myr 


2.28 


8 


6 


9 4K+0-03 
y - 4 °-0.31 


4 7n+ ' 19 

*• ,u -1.49 


8.96+2:* 




const 


3.54 


7 


8 


9 87 +o.io 
a -°'-o.io 


±z - oa -1.61 


8 QR+0.00 




const 


8.83 


9 


13 


q or+0.05 
aJO -0.10 


n nn+ ' 00 


7.911;?= 




burst 


6.01 


9 


IS 


8 99+° 22 
o.»»_q 24 


2 - 08 ^2'.08 


B.31+*;« 


r 


= 100 Myr 


1.43 


8 


2.3 


9 82+° 10 
a - oz -0.06 


o-oo±U% 


7 q fi +0'5° 
'' aD -0.05 




burst 


1.89 


9 


24 


10.25+°;°! 


00+° 00 
u - uu -o.oo 


a qi+0.00 
°- ol -0.00 




burst 


9.27 


9 


25 




00+ 7 77 
u - uu -o.oo 


7 qi+0.95 
'• al -0.05 




burst 


3.84 


9 


26 


9 34+0.11 


00+ ' 00 

uuu -0.00 


7 qg+0.10 
'■ au -0.10 




burst 


4.06 


7 


27 


9 20+ 020 

a - zu -0.09 




T RS+ ' 93 
<' bS -0.09 




burst 


1.25 


9 


157 


q 07+0. 16 
a - JI -0.20 




ok+0.20 
8 "™-g.40 




const 


10.65 


10 


248 


9 89+ - 08 

a - oa -0.14 


9.07±°J 


01 +0.05 
S ' 8i -0.10 


r 


= 500 Myr 


6.27 


10 



confidence limits marginalising over all other parameters if 
the degeneracies in the model set are dominating these un- 
certainties. For the errors in mass and SFR, a minimum 
error is determined from the range of normalisations of the 
best-fitting models that give Av; 2 <1. 

First, we see that with this restricted template set all 
of the objects are assigned acceptable fits at the 2a level 
(X 2 i$ 14 (12.5, 11 or 9.5) for 10 (9, 8, or 7) photometric 
points allowed a free choice of template, normalisation and 
age). 

Many of the objects are well fit by a single burst, in 
which case, if they have blue UV-slopes, the best-fit ages are 
quite low. This is because when objects with current star- 
formation are fit by a burst model, the age is determined by 
the rest-frame UV and, in the absence of star-formation in 
the templates, the UV must be fit by the young stars still 
present in a very young burst. The youngest stars die off 
over the first tens of Myr, building up a deep Balmer-break 
over a short amount of time. This allows the observed red 
H — 3.6/im colours to also be fit with young ages. Those 
objects that are best-fit by a burst but seem to show zero 
error in the SFR have no fits by constant or r models within 
A X 2 = 1. 

We see that the masses and SFRs are reasonably well 
defined, although some show improbably small errors, due to 
sparse sampling of the physical properties by this restricted 
model set. Masses are constrained to within a factor of ~ 2 
and SFRs seemingly better than this, although many of the 
best-fitting models are burst templates which, by definition, 
have no on-going star formation. The ages vary between 
~ 100 Myr to ~ 1 Gyr but the confidence intervals indi- 
cate that the ages are typically constrained to within a few 
hundred Myr. 

The results from this restricted model set therefore indi- 
cates that roughly a third of our sample harbour old stellar 
populations (> 500 Myr), and that the SFRs and masses 
appear fairly well constrained and are of the order 0-20 Mq 



6.2 Previous analysis (template set B) 

The ages derived in the previous section agree well with 
those reported by |Eyles et al.| |2007[ ) and four galaxies from 
of our sample were analysed in that work. Here we com- 
pare our results directly to their work, using our improved 
photometry but a similar template set. Table [4] shows the 
comparison between the results of our fitting analysis to 
theirs, restricting our template set to be representative of 
that used to derive the parameters published in Table 3 of 
Eyles et al. (2007); solar metallicity, no dust reddening, tau 
models with r < 1 Gyr, an instantaneous burst and a con- 
stant SFR model. |Eyles et &T\ ( |2007[ ) use a Salpeter IMF, 
so here we report their masses and SFRs re-scaled to the 
Chabrier IMF. 

The mass constraints for objects 6 and 23 agree very 
well. The masses for objects 5 and 8, which are the two 



objects analysed in Eyles et al. (20051, are both lower based 



on our fits, although still consistent within a factor of two. 

The differences between derived masses for any of the 
objects, but for objects 5 and 8 in particular, can be at- 
tributed to the differences in photometry. The CANDELS 
WFC3 data and deep A'— band data constrain the rest-frame 
UV to much higher accuracy. The deconfusion of the exist- 
ing IRAC data is also better constrained by using the high- 
resolution WFC3 image priors. 

The ages are mostly consistent, although object 6 is now 
fit with a considerably older template. The best-fitting ages 
of objects 6 and 8 are actually approaching the age of the 
Universe at z ~ 6 and would have been ruled out in the 



Eyles et al. ( 2007 I analysis 



The SED fits for objects 5 and 8 are plotted in 



These objects were first analysed in Eyles et al 



- 1 and M* = 10 9 - 10 10 M , respectively. 



who suggested that they harbour old stellar populations at 
high redshift (> 500 Myr). Using improved photometry but 
a similar template set we get consistent results, with the 
red rest-frame UV-to-optical colours being reproduced by a 
strong Balmer-break. The SFRs derived from both studies 
also show good agreement. 
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Table 4. A comparison of the results of Eyles et al. (2007) and best-fit parameters recovered here using a restricted template set (B), 
which was chosen to be consistent with the templates used in that work. Columns f-5 list the ID, best-fitting SFH, stellar mass, age and 
SFR based on SED fitting with template set B. Columns 6-fO list the values derived by Eyles et al. (2007) for the same objects. The 
best-fitting masses and SFRs derived in this work are corrected to total using the appropriate aperture correction. The reported values 
from Eyles et al. have been re-scaled to be consistent with a Chabrier IMF. 



ID 



model 



This work 
stellar mass 
(/10 10 M Q ) 



age 
(/Myr) 



SFR 

(/M e yr" 



1D 



model 



Eyles et al. 1 2007 1 

stellar mass age 
(/10 10 M Q ) (/Myr) 



SFR 

(/M Q yi- 1 ) 



5 
6 
8 
23 



const 
const 
= 1 Gyr 
= 100 Myr 



0.86 
0.31 
1.09 
1.03 



571 
905 
806 
286 



21.8 
5.2 

13.3 
9.0 



31.2185 
23_2897 
23.6714 
32.4331 



const 
= 1 Gyr 
= 500 Myr 
= 70 Myr 



1.3 
0.4 
1.7 
0.9 



640 
640 
720 
260 



23.8 
5.5 

13.3 
4.4 



A (//im) 



0.6 



0.4 



0.2 



0.0 



0.0 \|/ 



- 




irr 


1 






1 ID 5 571 Myr 


E(B-V)=0.0 


in Jm 









0.2 



0.4 
A, (/Vm) 



0.6 



0.E 




0.4 

K (/^ m ) 



0.8 



Figure 2. The best-fitting SEDs to objects 5 and 8, with a template set restricted to solar metallicity, zero reddening models with 
r ^ 1 Gyr, constant SFR or an instantaneous burst. These two objects were first analysed in Eyles et al. (2005) who suggested that 
they harbour old stellar populations at high redshift, with their analysis suggesting ages > 500 Myr. Using improved photometry, but a 
similarly-restricted set of templates, our results agree. 



Table 5. Results of SED fitting with the full set of smoothly- varying SFHs, extinction and metallicitics (template set C). The first column 
gives the IDs of the objects. Columns 2-7 give the best-fitting stellar mass, SFR, sSFR, extinction, metallicity and age respectively. The 
sSFR for burst models (SFR = Meyr -1 by definition) are estimated using the SFR of the next best-fitting model that is not a burst. 
All quoted errors are the 68% confidence limits, marginalising over all other free parameters. Column 8 gives the lower limit to the Age 
from the la uncertainties from all the models except the burst. The final two columns give the best-fitting SFH and \ 2 - The typical 
degrees of freedom for 9 photometric points is 4. 



Object 


logirj(Mass) 


SFR 


sSFR 




E(B-V) 


metallicity 


logio(Age) 


Age lower 




Model 


x 2 


ID 


(Ms) 


(/Moyr- 1 ) 


(/Gyr- 1 ) 






(/z©) 


(/Gyr) 


limit (lcr /Myr) 








2 


9.50t°,; 22 


00 + 27 ' 6 
u ' u -0.0 


1.2 





uuo -0.005 


0.4 




32 




burst 


2.79 


5 




19 6 + 40 2 


2.8 





uuu -0.000 


0.4 




43 


EI 


sSFR.=2 Gyr- 1 


2.11 


6 




4.7 til 


1.7 





mo +0.044 
uuu -0.000 


0.2 


8.96+^5 


128 




const 


3.54 


8 


n nc+0.12 
a ' MD -0.14 


199 +3.9 
-2.6 


1.3 





m9 +0.032 
uuz -0.000 


0.2 


8.96+°'°° 


500 




r=2 Gyr 


7.06 


13 


q 14 +0.45 


n+ 132 ' 5 

u ' u -0.0 


2.3 





1K „+0.136 
loo -0.091 


1.0 


v 19 +l!84 
' - iz -0.10 


13 




burst 


2.82 


IS 


a ' U -0.54 


+ 45 ' 7 
uu -0.0 


2.7 





m[ -+0.279 
uuo -0.005 


0.02 


7 SK+ 1 ' 10 
'■ St> -0.84 


10 




burst 


1.23 


23 


9.90l° ;| 


n o + 96 - 5 

u ' u -0.0 


0.3 





nl0 +0.210 


0.02 


8.06tg;gj 


32 




burst 


1.75 


24 


io.o6±g:« 


n n +506.5 
u ' u -0.0 


2.3 





qci +0.020 
ool -0.175 


0.2 


7 02+ 1 ' 99 
' ■"''-o.oo 


10 




burst 


1.23 


25 






1.7 





n42 +0.153 
u * z -0.042 


1.0 


„ 71 +0.30 
S -' i -1.24 


20 




r=0.5 Gyr 


2.28 


26 


9 18+ ' 52 


+137.4 
u ' u -0.0 


2.4 





1sn +0.l60 


0.4 


7 99 +1.74 
1 ■''''-0.20 


10 




burst 


1.93 


27 




3.5 + 3 H 


1.6 





mo +0. 151 
uuu -0.000 


0.4 


8.l6tg;g 


10 




r=0.05 Gyr 


1.02 


157 


9.43±°;f 5 


i5.o + 7 -j 


5.6 





002 +0.064 
uuz -0.002 


0.2 




32 


EI 


sSFR=4 Gyr- 1 


8.03 


248 


9.89±g:i§ 


9.6 i 3 3 6 


1.2 





nm +0.030 
uuu -0.000 


0.2 


8.91+SJS 


363 




r =1 Gyr 


6.23 
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6.3 All smoothly varying SFHs (template set C) 

Although all of the objects are assigned acceptable fits from 
the restricted template set, we now look at how the de- 
rived properties and their constraints change when we allow 
extinction by dust, multiple metallicities and an expanded 
range of SFHs (including exponentially increasing SFR mod- 
els). 

The best-fitting parameters and 68% confidence limits 
are reported in Table [H] In this case, we require x 2 J; 11 
(9.5, 7.8 or 6) for the fit to be considered acceptable with 10 
(9, 8, or 7) photometric points. Again, all objects have been 
assigned an acceptable fit and the x 2 values are considerably 
improved for most objects. In allowing for dust in the fitting 
to these objects, the uncertainties in SFR and age increase, 
as expected due to the known degeneracies in these param- 
eters. The masses derived from the restricted model set are, 
however, proven to be quite robust. 

In general the masses, ages, extinction and metallicity 
tend to agree between the different models, although burst 
models tend to give consistently lower best-fitting masses 
and ages, primarily because younger ages are needed to fit to 
the rest-frame UV in the absence of on-going star formation 
in the templates, as described in section 6.1. We can see this 
agreement when plotting Ay 2 surfaces for various parameter 
combinations for a representative object (object 2) in Fig [3] 
The la limit is plotted in pink and shows the 68% joint 
confidence interval (Ay 2 < 2.3). These plots show that the 
parameter ranges overlap very consistently for the different 
models. 

The differences between the limits can be described 
mainly by the differences in the allowed ranges of param- 
eter space occupied by different models. For example, mod- 
els with a constant star-formation rate must have a mini- 
mum normalisation required to match to the luminosity of 
the sources, thereby the lower SFRs permitted by r and EI 
models are not seen here. The agreement in overall normal- 
isation of the models is what drives the overall agreement 
in derived masses. The resultant SFRs are therefore driven 
by the range of SFRs permitted within this normalisation 
which can also reproduce the galaxy's observed rest-frame 
UV. 

Ages are not well constrained, as expected, because ei- 
ther age or dust reddening can produce the same observed 
SED, for any SFH that is not increasing with age. The de- 
generacies are slightly different for the EI models, due to 
older galaxies having higher star- format ion rates and we see 
no apparent correlation between age and reddening in the 
likelihood contours (Fig. [3]F) - 

Most of the galaxies allow for dust within the confi- 
dence contours, whether the best-fitting model requires it 
or not. |Eyles et al.| ( |2007[ ) take the absence of dust from 
the best-fitting models to show little evidence for intrinsic 
reddening within the galaxies, although their derived ages 
range from ~ 300 - 700 Myr, requiring older, redder popu- 
lations to fit to the observed SEDs. As we have seen from 
restricting our model set, our results agree. When looking at 
all the models, metallicities and reddening values, however, 
a wide range of ages provide acceptable fits to each object 



(also remarked on by |Eyles et al. ( 2007 1 with reference to 
the minimum and maximum ages derived from burst and 
constant SFR models), and some objects do show evidence 



for intrinsic reddening in their SEDs (see Fig. [5]l . We also 
note that the older age estimates are at odds with the predic- 
tions of various models (e.g. Finlator et al. ( |2010| ) predict 
ages ranging from a few Myr to ~150 Myr at these red- 
shifts). Consequently it is of interest to investigate whether 
any of the objects actually require an old stellar population 
(> 300 Myr) which dominates the stellar mass. In Section 
6.3.2, we consider again the degeneracies inherent to our lim- 
ited model set and consider the full range of possible models 
within the marginalised 68% confidence contours in order to 
see whether these old ages are truly robust. 

6.3.1 Exponentially increasing SFHs 

As described in Section 4.1, the exponentially increasing 
SFHs are initiated with a seed mass produced by a burst 
of initial star formation. This is just a by-product of the 
parameterisation, but this seed-mass leaves an imprint on 
the observed SED by reddening the observed UV-to-optical 
colours. 

We see that two of the objects are best-fit by EI SFHs, 
although for a constant sSFR of 2 Gyr -1 a galaxy needs to 
have been growing in mass over the last ~ 890 Myrs to make 
the contribution to the SED of the original burst insignif- 
icant (i.e. < 10% rest-frame optical flux in the absence of 
nebular emission). As a consequence, the best-fitting model 
for one of these objects is essentially an instantaneous burst 
with some ongoing star-formation. 

6.3.2 Template degeneracies and age determination 

Our results with smoothly-varying SFH models show that 
three of the sample of 13 galaxies have old best-fitting ages 
(> 300 Myr, objects 6, 8 and 248) and all of the objects 
that are not best-fit by a pure burst template show best- 
fitting ages of at least 200 Myr. As discussed previously, the 
burst templates give very young best-fit ages when fitting 
to objects that do show star formation, simply because the 
youngest bursts still have young, hot stars that can fit the 
rest-frame UV flux. We therefore look to the lower limits in 
age (lowest ages allowed within Ay 2 = 1), when burst mod- 
els are not included in the template set, which are reported 
in column 8 of Table[5] We see that, when the SEDs are well 
described by a burst, they also allow for very young ages 
when the burst models are not present, because very young 
exponentially increasing SFHs also show a strong burst com- 
ponent, as discussed above. 

As for the objects with best fits which indicate an older 
population, the lower limits in age show that two of the 
objects, 8, and 248 appear to require a population of at 
least ~300 Myr in their fits. 

Within this context it is interesting to examine the fits 
for typical examples of objects with consistently old and 
young best-fits (Figs [4] and [5] respectively) . For most of the 
SFHs (except EI), the acceptable models with dust have 
the youngest ages (Fig. [5| while the objects being fit with 
consistently-old stellar populations (> 500 Myr) have bluer 
UV slopes that allow for less intrinsic reddening (Fig. Ek. 
While the objects with the youngest best-fitting ages may 
still have fairly red UV-to-optical colours, they tend to have 
redder UV-slopes and evidence of intrinsic reddening from 
their SEDs. 
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Figure 3. Fitting to object 2: log(stellar mass) vs log(SFR) of the models (a-c) and model age vs. extinction (d-f) for three different 
sets of models; r models (exponentially decreasing SFR), constant SFR and exponentially increasing SFR. The greyed areas show the 
parameter space covered by the best-fitting models. In the case of age and extinction, these are supplied as a grid of values, while both 
SFR and mass are the values of the best-fitting models within the grids of parameters. The pink, green and blue regions show the ranges 
of parameters with best-fitting model \ 2 values within A\ 2 = 2.3, 6.2 and 11.8 respectively, of the overall minimum x 2 value. 



Fig. [6] (a) and (b) show the UV-to-optical (H — 3.6^m) 
and UV (J—H) colours of the objects plotted against the age 
of the best-fitting model. Again, we see that the objects that 
are assigned the oldest ages from this template set tend to 
have fairly red UV-to-optical colours, as would be expected 
if the optical fluxes are dominated by older stellar popu- 
lations, but they also have very blue UV slopes with very 
little dust allowed in the fitting (filled symbols). We note 
that the object with the youngest best-fit age (object 24) 
also has a very red H — 3.6/^m colour, yet an older stellar 
population would not have time to build up for such a young 
galaxy. This object also has a red UV-slope, however, and 
the template requires a relatively high level of extinction 
(E(B— V)=0.351, see Table 5) which explains how an object 
with such a red H — 3.6/im colour can be fit with such a 
young model. 

Those SEDs that allow for the presence of dust in the 
shape of their UV-continuum therefore allow for a wider 
range of ages in the fitting, generally extending the allowed 
range to younger models, whereas a blue UV-continuum plus 
a red UV-to-optical colour allows for less freedom when ap- 
plying dust. The resultant best fit is then a balance between 
the age required to fit the UV-optical plus current star for- 
mation to fit to the UV. 

These slightly counter-intuitive results, with older ages 
assigned to objects with evidence of current star-formation, 
are indicative of a conflict between the observed SEDs and 
the model templates when blue UV slopes are combined with 
apparently large Balmer breaks. This is reminiscent of the 



difficulties reported by Labbe et al. ( 2010 \ in simultaneously 



fitting to a blue UV-slope and red UV-to-optical colours. In 
Section 7 we proceed to investigate whether a two compo- 
nent SFH and the addition of nebular emission can resolve 
this tension. 



7 SED FITTING - TWO COMPONENT STAR 
FORMATION HISTORIES (TEMPLATE 
SETS D-E) 

In this section we explore the possibility of explaining the ob- 
served SED characteristics with a more general parameteri- 
sation of the SFH and, in the process, see whether this solves 
the conflict between blue UV-slopes and red UV-to-optical 
colours. This parameterisation uses a two-component model 
that decouples the instantaneous SFR from the mass built 
up over time. The results both with and without nebular 
emission are studied (template sets E and D respectively). 
These models consist of a constant duration of star forma- 
tion to describe the current star-forming component and an 
older burst to fit the rest-frame optical (a full description is 
provided in Section 4). 

The recent durations of star formation are fairly 
sparsely sampled with only three different durations (10, 30 
and 100 Myr), but by parameterising the star- forming com- 
ponent in this way, rather than by a second burst, we can 
derive the SFRs more consistently from the models, with- 
out using calibrations between UV and SFR based on spe- 
cific SFHs that do not hold for a very recent burst of star 
Madau et al.|1998[ ). 



formation (i.e. 



Ideally, this parameterisation will give upper and lower 
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Figure 4. SEDs and best-fitting models of three galaxies that 
apparently require old stellar population ages to provide good 
fits. The red points plus error bars show the observed SED and 
the blue points show the synthetic photometry based on the best- 
fitting model. 



limits to the mass (Papovich et al. 2001 1 and age of the 



dominant stellar component (i.e. the age of the burst with 
minimum of 85% of the mass according to this parameterisa- 
tion) . We have already seen that pure bursts provide a lower 
limit in age when the UV indicates the presence of young 
stars. To be able to place improved constraints on the instan- 
taneous SFR, however, an extra constraint is added to the 
fits. The models are also required to produce at least enough 
ionising photons to explain the Lya emission for all objects 
except the two for which we were unable tie the calibration 
of the spectroscopy to the photometry (objects 18 and 16, 
see Section 2.3). This constraint ensures that a very young 



2 K (/m™) 




0.4 

K (/^ m ) 

Figure 5. As for Figure |4]but for objects which are consistently 
best-fit by young (< 300 Myr) stellar populations. The grey point 
in the plot for Object 26 shows the filter excluded from fitting due 
to the uncertainty in the Lya contribution. 



burst cannot dominate the UV part of the spectrum when 
the observed Lya indicates that there should be a contribu- 
tion from current star formation. This is important for the 
SFR and age estimates as it prevents the lower limits being 
set by models in which the UV flux is dominated by young 
stars from a young burst, rather than from the star-forming 
component. 

A lower limit on the flux of ionising photons is calcu- 
lated using the calibrated Lya flux, which is in any case a 
lower limit of the intrinsic Lya due to absorption by the 
IGM, and assuming Case B recombination | |Osterbrock fc| 
Ferland|2006 1 giving: 

(1) 



L(Lya) = 24 x 4.87 x 10~ 13 (1 - f BS c)N L yC 
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Figure 6. The best-fitting age vs H — 3.6/im colour, which spans 
the 4000A/Balmer break region (a), and the J — H colour, which 
samples the near-UV slope (rest-frame) (b). Filled symbols show 
best fitting solutions which only allow very small amounts of dust 
(E(B— V) < 0.05 within la uncertainties). The three objects with 
the reddest J — H colours also have the youngest best-fitting ages 
as they have been fit with burst templates. The observed correla- 
tion between J — H colour (rest-frame UV-slope) and model age 
may be an indication of inherent model degeneracies producing 
older age estimates in objects that allow for little reddening by 
dust. 



where Nl V c is the flux of Lyman continuum photons, 
L(Lya) is the luminosity of Lya in erg s and f esc is set 
to to provide the lower limit. 



7.1 Templates without nebular emission 
(template set D) 

7.1.1 Results 

Table [6] documents the best-fitting parameters with the 
la limits marginalised over all the different models within 
the set. In general the derived masses and SFRs agree 
well with the confidence contours derived from the full set 
of smoothly-varying SFHs, although when the objects are 
best described by a burst in template set C, their best- 
fitting masses are in general lower than when fit with the 
burst+constant models here. 

Those SEDs with redder UV colours allow for more ex- 



tinction and so the uncertainties in derived SFR are larger, 
but this is also true when fitting them with the smoothly 
varying SFHs. We find that two of these objects that allow 
for significant extinction (ID 13, 24) have very high best- 
fitting SFRs (83 M Q yr" 1 and 438 M yr" 1 respectively). 
These high SFRs are due to the rest-frame UV and most 
of the rest-frame optical flux being dominated by the star- 
forming component, with the older burst component work- 
ing to boost the flux in the optical. However, it is worth 
noting that younger burst components in combination with 
lower constant SFR components are also allowable within 
the uncertainties. 

Many of the best-fits indicate old burst components 
with small star-forming fractions (< 3%). Firstly, taking 
this at face value, we can say it's incredibly unlikely that 
these objects formed most of their stars in a large starburst 
of short duration at very high redshift, given current ob- 
servational constraints. For example, a burst just 0.1 Myr 
old, with 0.02 Z metallicity, that produced 10 9 M of 
stars would be brighter than Hab = 24.4 when observed 
at z ~ 8, although extinction with E(B — V) = 0.35 could 
bring the brightness below the detection threshold for the 
CANDELS deep //—band imaging (5a point source depth of 
27.7, |Grogin et aL]|2011[ ). It is possible, however, that these 
objects could be described by stochastic star formation his- 
tories that would leave an imprint on the observed SEDs 
that would be similar to an old, instantaneous burst. 



7.1.2 New age determinations and intrinsic uncertainties 

In this section we investigate the ages derived with the 
burst+constant models to see whether any of the objects 
still require an old stellar population once the instantaneous 
SFR has been decoupled from the accrued mass. As we have 
set the constraint that the model must supply enough ion- 
ising flux to describe the observed Lya, the lower limits in 
age (the youngest models supplying fits within Ax 2 = 1) 
are not just set by the age of a burst that fits the rest-frame 
UV, but must also supply the minimum flux of ionising pho- 
tons required to explain the Lya emission. However, caution 
must be used regarding the results for objects 18 and 26 be- 
cause we were not able to apply this constraint and we see 
that the lower age limit is set by the age of a pure burst, 
with no current star-formation. We know these two objects 
have Lya in emission so these lower limits are likely to be 
unphysical. 

Many of the objects do not have an upper age con- 
straint, with the upper limits butting up against the age 
of the Universe. This is due to the degeneracy between the 
star-forming fraction and the age of the burst (as well as, 
to a lesser extent, extinction and metallicity). Essentially, a 
model with a fairly young burst (> 100 Myr), which consti- 
tutes a majority of the flux in the rest-frame optical, plus a 
certain star-forming fraction, looks very similar to a model 
in which the star-forming fraction dominates the rest-frame 
optical flux but an old, massive burst also makes a signifi- 
cant contribution. In other words, the older a population is, 
the easier it is masked by any current star-formation. We 
cannot, therefore, rule out old ages for most of the objects 
although we provide the caveat that the older the burst, the 
larger the mass of stars required to be built in a shorter 
amount of time, and hence the less likely it is to represent 
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Table 6. Results of SED fitting with the set of burst+constant star formation models (template set D). The first column gives the IDs 
of the objects and the second indicates the duration of constant star formation of the best-fitting model in Myr. Columns 3-10 give the 
best-fitting mass, SFR, sSFR, extinction, age, la lower limit in the age, fraction of mass in the star-forming population and metallicity 
(as a fraction of Solar) respectively. All quoted errors are the 68% confidence limits, maginalising over all other free parameters. The final 
column gives the x 2 of the best-fitting model. The typical degrees of freedom for for 9 photometric points is three. 
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reality. Instead we continue to ask the question of whether 
old ages are required to fit to the observed SEDs. 

First we see that some of the objects do not have any 
true age constraints, with lower limits defined by the lowest 
ages of the models and upper limits approaching the age 
of the Universe (objects 13, 18, 24 and 26). These either 
have no constraint set by Lya, either with poorly calibrated 
Lya emission (26) or because the system is an absorption 
system (24), or the SEDs allow for a lot of intrinsic reddening 
(13, 24, 26). We cannot constrain the age when reddening 
introduces so many degeneracies or if fits are allowed the 
freedom to fit the majority of the UV with the non-star 
forming burst component. 

The two objects that required an old population from 
template set C, objects 8 and 248, have constrained lower 
age-limits of 115 and 290 Myr respectively. Are these ages 
robust when nebular emission is added to the templates? 



7.2 Templates with nebular emission (template 
set E) 

The template set with full nebular emission is produced 
by adding nebular emission (continuum and line emission) 
to the burst+constant star formation templates, using the 



method of Robertson et al. (2010) which is summarised in 



McLure et al. (2011). As we are primarily interested in the 



effect of nebular emission on the derived ages and masses 
of the objects, we restrict ourselves to a single value of 
the escape fraction of ionising photons, f esc — 0.2, which 
is motivated by the typical value of / esc which is required 
for star-forming galaxies to maintain reionisation at z ~ 7 



(Robertson et al. 20101 



Nebular continuum and line emission is then added to 
the burst+constant templates using the number of ionising 
photons per second, Nl V c, provided as an output from the 
BC03 models. This number is used to calculate the H/3 lu- 
minosity, and other HI line intensities follow from the ratios 



predicted by standard recombination theory ( Osterbrock & 



Ferland 20061. Emission lines of metals are added using 



the empirically- derived, metallicity-dependent, line intensi- 



ties (relative to H/3) from Anders & Alvensleben (20031 



7.2.1 Results 

We find that a subset of the objects (IDs 5, 6, 8, 18, 27, 
157, 248) give marginally better fits when nebular emission 
is added to the models, although for most of these objects 
the x 2 improvement is not statistically significant (the min- 
imum \ 2 values are within A\ 2 = 1). We cannot say that 
these objects require nebular emission from the templates to 
provide suitable fits, but we can investigate how the inferred 
physical properties are affected as well as the properties of 
the SEDs that are better reproduced when nebular emission 
is included. We defer discussion of the second point to Sec- 
tion 8.1 and concentrate here on how the inferred physical 
properties of the objects can change. 

The best-fitting models from this template set give sys- 
tematically lower best-fitting masses, with a median offset of 
0.18 dex, compared to the same template set without neb- 
ular emission. This is comparable to the results of |McLure] 
et al. (20111 where, for a subset of their photometrically- 
selected 6 < z p hot < 8.7 LBGs detected at 3.6/im, 16 out of 
21 objects show an offset of less than 0.4 dex when compar- 
ing results with and without nebular emission (the templates 
used for fitting to the object SEDs consisted of a single ex- 
ponentially decaying SFH with two possible metallicities). 
The lower mass estimates are caused by additional nebular 
line contribution to the optical fluxes. A burst with a deep 
Balmer break is, in many cases, no longer needed to fit to 
the observed optical fluxes, as found by other studies (e.g. 
Schaerer fc de Barros|2009| |Labbe et al.||2010 1. 

The best-fitting SFRs derived with nebular emission, 
however, are very similar to those derived from the tem- 
plates without nebular emission. This is due to the rest- 
frame UV constraining the amount of star formation and 
the inclusion of nebular continuum emission, which acts to 
redden the rest- frame UV, only minimally affecting the SFR 
determination. 

The SFRs are actually better constrained for this tem- 
plate set compared to sets D (and C) because the higher 
mass, higher SFR models no-longer provide adequate fits to 
the SEDs. This is because a high SFR produces a large con- 
tribution to the rest-frame optical fluxes via nebular emis- 
sion lines that were not accounted for in the previous tem- 
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plate sets. Nebular emission therefore limits the contribu- 
tion of the burst component to the IRAC fluxes, therefore 
limiting the derived masses and SFRs. For example, the ex- 
tremely high SFRs from template set D are now not allowed 
within A% 2 = 1 for objects 13 and 24 due to the addition 
of nebular emission, although the best-fitting SFR of object 
24 is still quite high. It is worth noting here that we assume 
the old and young stellar populations show the same levels 
of extinction and that the nebular extinction is the same 
as the stellar extinction. Relaxing these constraints would 
likely show that the SFRs were as uncertain as when esti- 
mated from templates without nebular emission, as it would 
break the direct correspondence between SFR and nebular 
line contribution to the optical fluxes. 

7.2.2 Age determination from templates with Nebular 
Emission 

Fig.[7|shows the best-fitting SEDs for each of the objects de- 
rived from each of the template sets, C, D and E. In particu- 
lar, the lower right-hand panel shows the best-fits to object 
248 which previously appeared to require an old stellar pop- 
ulation. The lower age limit for 248 is still fairly high at 227 
Myr, although it can be seen that its blue 3.6/im— 4.5/im 
colour is better fit by the models including nebular emis- 
sion, they have some difficulty in reproducing the very blue 
colours. This is discussed in more detail in Section 8.2. 



8 PROPERTIES OF UV-BRIGHT LBGS 

8.1 Observed SED characteristics driving inferred 
physical parameters 

We have used a number of different template sets for SED- 
fitting to our sample of spectroscopically confirmed LBGs. 
The differences between the inferred physical properties 
from each template set illustrates how restraining our priors 
in SFHs, extinction, metallicity and nebular emission can 
influence our conclusions. As has been well established, the 
most uncertain parameters are the ages and metallicities of 
the objects. When nebular line emission, in particular, is 
omitted from templates then the inability to pin down the 
age and metallicity of the object does not impact signifi- 
cantly on the derived masses and allowing for a full range of 
each in the fitting allows for a realistic determination of the 
uncertainties. 

When nebular line emission is added, however, we see 
that it is possible for a recent starburst with strong nebular 
emission lines to mimic the signal of an old stellar popula- 
tion. For certain redshift and filter combinations (in partic- 
ular, for our sample at z ~ 6 with the only two filters deep 
enough to sample the rest-frame optical region of the SED 
being the IRAC 3.6 and 4.5/tm channels) this can lead to 
a case where the nebular emission contribution cannot be 
constrained, or equivalently the contribution of any under- 
lying old population is also not well constrained, leading to 
marginally larger uncertainties on the derived stellar masses. 

The SFRs derived from different template sets show 
good agreement, in general, because the deep near-IR imag- 
ing provides strong constraints in the rest-frame UV. For 
objects that show evidence of intrinsic reddening in their 



SEDs, the SFRs are unconstrained because the level of ex- 
tinction is very uncertain. Adding nebular emission to the 
templates limits the maximum estimates of SFR as they 
would be required to contribute to the rest-frame optical 
fluxes via nebular line emission, but this relies on our as- 
sumption that the nebular and stellar components suffer the 
same extinction. 

The observed SEDs and best-fitting model SEDs for the 
template sets C (all smoothly varying SFHs), D and E (two 
component SFHs without and with nebular emission, respec- 
tively) are plotted in Fig. [7] We can see that there is very 
little difference between the fits derived from any of the tem- 
plate sets, with the small errors from the deep HST imaging 
providing very good constraints in the UV. There are differ- 
ences in the contributions to the optical fluxes, although for 
the most part they are also marginal and in general can be 
attributed to the poorer constraints. For the objects show- 
ing blue 3.6/mi — 4.6/im fluxes, however, this is not the case 
and they are slightly better reproduced by the templates 
including nebular emission. This may indicate that there 
is indeed nebular emission contributing to the broad-band 
fluxes of these objects but the improvement in the fits are 
not statistically significant. 

Plotting the colour-colour space covered by the models 
and observations for a range of different SFHs in Figs [8] and 
[9] we see that templates including nebular emission (dark 
grey regions) extend to bluer colours in 3.6/im — 4.5/im 
because, at this redshift, [OKI] and H/3 are contributing to 
the 3.6/im flux and Ha is contributing to the 4.5/im flux. 
None of the models are able to reproduce blue IRAC colours 
with the appropriate range in H — 3.6/im colour even with 
a contribution from nebular emission lines. Even allowing 
for ages <10 Myrs in the smoothly varying SFHs, the bluest 
observed 3.6 — 4.5/im colours are not reproduced by any of 
the SFHs plotted in the figure. A potential reason for this is 
discussed further in the next section. 

8.1.1 IRAC colours 

Fig. [8] (c) shows that the two-component models follow the 
observed correlation between the 3.6/im — 4.5/im and J — H 
colours. This correlation is expected when considering red 
3.6/tm — 4.5/tm and J — H colours, as extinction would act 
to redden both the rest-frame UV and the rest-frame optical. 
In this projection, we can see that the models both with and 
without nebular emission do quite well at reproducing the 
blue 3.6/im — 4.5/im colours. 

Fig.|8|a) shows the 3.6/tm— 4.5/im vs. H— 3.6/im colours 
for the two component models both with and without neb- 
ular emission (dark grey and blue regions respectively) with 
the observed colours over-plotted. The observed 3.6/im — 
4.5/tm colours extend to much bluer values than covered by 
the models. Although it is possible that there may be sys- 
tematic errors introduced by the IRAC deconfusion that are 
not included in the errors plotted, it is important to question 
whether something is missing from the models as the high 
3.6/im flux could be compensated for by allowing a larger 
Balmer-break from an old underlying stellar population, as 
we found to be the case when nebular emission is not in- 
cluded in the templates. 

One plausible limitation of the nebular models is that 
the average line ratios used to assign line fluxes to colli- 
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Table 7. Results of SED fitting using the set of burst+constant star formation models with added nebular emission (template set E). The 
first column gives the IDs of the objects and the second indicates the duration of constant star formation of the best-fitting model in Myr. 
Columns 3-10 give the best-fitting mass, SFR, sSFR extinction, age, lower age limit, fraction of mass in the star-forming population and 
metallicity (as a fraction of Solar) respectively. All quoted errors are the 68% confidence limits, maginalising over all other free parameters. 
The final column give the best-fitting model and \ 2 - The typical degrees of freedom for 9 photometric points is 3. 
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sionally excited lines in the nebular templates may not be 
representative of the ratios expected for these objects. In 
particular, a higher f([OIII])/f(H/3) ratio would boost the 
[OIII] flux relative to Ha and hence produce bluer IRAC 
colours. Rest-frame optical spectroscopy of LBGs at lower 
redshift (z ~ 3) give a spread in measured f([OIII])/f(H(]) 
ratio, although the derived metallicities are all close to so- 
lar (Mannucci et al. 20091. Finkelstein et al.| ( |2011[ > also 
obtained spectroscopy of two LAEs at z ~ 2.3 that show 
strong [OIII] emission, with no indication of any AGN con- 
tribution to the ionising flux. Given the strong temperature 
dependence of [OIII] emission compared to that of H/3, this 
ratio would be sensitive to HII region geometry within the 
galaxies and so could present a wide range of values within 
the population. 

In Fig. |8ja) we illustrate how the IRAC colours are 
affected by boosting the [OIII] flux of the 0.2Zq metallic- 
ity models (the metallicity that gives the highest intrinsic 
f([OIII])/f(Hfl) ratio) by a factor of 1.5 giving a flux ra- 
tio comparable to the highest flux ratio measured for the 
objects at z ~ 3 (Mannucci et al. 2009 1 . This is purely 



5.5 x 10~ 17 erg s _1 cm~ 2 (Xue et al. 20111 and, following 



the argument laid out in Vanzella et al. (20101, this can be 



illustrative, as without rest-frame optical spectroscopy we 
cannot constrain the f ([OIII])/ f (H/3) ratio at these red- 
shifts, although it can be seen that boosting the [OIII] flux 
in this way does improve the coverage of colour space of the 
models. 



8.1.2 X-ray emission 

Given the possibility that the blue observed 3.6/^m — 4.5/im 
fluxes can potentially be explained by boosted [OIII] EWs, 
it is important to consider the likelihood that any of these 
objects harbour an AGN that could be driving the strong 
[OIII] emission. First we note that there has, as yet, been 
no direct evidence of an AGN contribution from spectro- 
scopic observations of low-redshift LBGs that also show high 
f([OIII])/f(H/3) ratios, as discussed in the previous sec- 
tion. 

We further note that none of these objects have matches 
in the 4Ms Chandra point source catalogues. The new Chan- 
dra 2—8 keV X-ray data has on-axis, S/N>3 sensitivity of 



used as an upper limit to the AGN-driven Ha flux using 
the observed correlation between X-ray luminosity and Ha 
luminosity ( Panessa et al.|2006 1. 

The limiting Ha luminosity can be related to 
the observed 3.6/im — 4.5/^m colours by finding the 
f ([OIII])/ f (H/3) ratio that would be required to boost the 
flux at 3.6/xm, assuming the continuum is flat and the lines 
are not affected by extinction (extinction would only act 
to redden the derived 3.6//m — 4.5/xm colours). We assume 
standard line ratios of f(Ha)/f(H/3) = 2.85 (case B recom- 
bination, ?) and /([O7JI]A49 59A)//([O77I]A5007A) = 2.9 
( jAnders fc Alvensleben|2003[ close to the measured value of 
2.99 for a sample of 62 AGN by |Dimitrijevic et al.|2007 l. 

The correlation between Ha flux and X-ray flux indi- 
cates that if the objects all had X-ray fluxes at the sensi- 
tivity limit, their AGN-driven Ha fluxe^] would be of order 
2.3 x 10 18 erg s" 1 cm' 2 . The f ([OIII])/ f (H/3) ratio that 
would be required to reproduce the blue 3.6/im — 4.5^im 
colours for four of the objects (IDs 2, 18, 157 and 248) are 
shown in Table[8] where the errors in / ([OIII])/ f (H/3) ratio 
are derived from the measured photometric errors. 

The required f ([OIII])/ f (H/3) for objects 2, 157 and 
248 are, in particular, high compared with measurements 
of the order ~ 10—15 for local AGN (Kauffmann et al 



|2003[ ). Although it is possible that a central AGN may be 
obscured in the X-ray, allowing for higher AGN-driven line 
fluxes for the Chandra 4Ms flux limit and hence a lower 
f([OIII])/f(H/3) ratio to produce the observed 3.6/^m — 
4.5/xm colour, this would still require the objects to have X- 
ray fluxes sitting at the flux limit of the Chandra survey. It 
is therefore unlikely, if the blue 3.6^tm — 4.5^tm colours are 
driven by emission lines, that the dominant ionising flux is 
supplied by a central AGN. 



3 The derived correlation between X-ray luminosity and Ha lu- 
minosity used here is taken from the total Panessa et al. (2006) 
sample plus low-redshift bright Seyfert 1 galaxies (Tot+QSO sam- 
ple, see their Table 3). 
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Figure 7. The best-fitting SEDs for each object from the following template sets: Smoothly varying SFHs (template set C, plotted as 
red diamonds), two component SFHs without nebular emission (template set D, plotted as blue triangles) and with nebular emission 
(template set E, plotted as yellow squares). The observed SEDs and corresponding errors are plotted as black circles, with any filters 
that have been excluded from the fitting due to an uncertain contribution from Lya being plotted in grey. 



8.2 Ages 

The median best-fitting ages derived for each of the template 
sets are 145 Myr (template set C), 182 Myr (D) and 288 Myr 
(E) respectively. The median age derived using template set 
C is lower than for the other two template sets because just 



under half of the sample are best-fit by young bursts. Fig. |10| 
shows the constraints placed on the ages for each object 
from the two-component templates combining the results 
obtained with and without nebular emission. 

The three objects with the best age constraints are ob- 
jects 2, 5, and 23 which show peaks covering the age range 
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(a) 



(b) 



(c) 



Figure 8. Observed colours of our z ~ 6 LBG sample plotted against the two-component model colours. From left to right the plots 
show: The IRAC 3.6/im - 4.5/^m vs. H — 3.6fim colours, the J — H vs. H — 3.6/^m colours and the 3.6/^m — 4.5/im vs. J — H colours. 
Objects that are better fit by models without nebular emission are plotted as red triangles and the objects that are better fit by models 
with nebular emission are plotted as blue diamonds. The filled symbols represent objects which were originally confused in the IRAC 
images, whereas the open symbols show objects that are relatively isolated. The dark grey regions show the colours covered by the models 
with nebular emission and the light blue regions show the colours covered by the models without nebular emission. The model colours 
are all plotted at z = 6, although they are minimally affected by allowing the full range of redshifts within our sample (5.5 < z < 6.2). 
The lighter grey regions show how the colours are affected by boosting the [OIIIJ/H/3 ratio (see text for details). 




(a) Constant SFR models (b) EI models (c) r models 




0.0 0.5 1.0 1.5 2.C 0.0 0.5 1.0 1.5 2.C 0.0 0.5 1.0 1.5 2.C 

H - H - 3.6/im H - 3.6/zm 



(d) Constant SFR models (e) EI models (f) r models 

Figure 9. The 3.6^tm — 4.5/^m vs. H — 3.6/^m and J — H vs. H — 3.6/xm colours of the models with smoothly-varying SFHs at z = 6, 
both with (grey regions) and without (blue regions) nebular emission. The darker blue regions show where these colours overlap. The 
observed colours are plotted as described in Fig. [8] The plotted colours include very young (< 10 Myr) populations. 



~ 50 — 200 Myr. The age constraints are, however, quite 
poor for many of the objects in the sample, in particular 
those objects without a full complement of filters (6, 18, 
26), which have insufficient data points once these filters are 
removed from the fitting to provide meaningful constraints 
(number of filters ^ 7). Without the 4.5/im filter, the neb- 
ular emission for object 6 is completely unconstrained, and 
without the filters that have contributions from Lya, as well 
as a minimum SFR constraint, the other three objects have 
very flat probability density functions (PDFs). 

The old median ages (~200 Myr) derived from the dif- 
ferent template sets are at odds with the extremely young 
ages often found from fitting with smoothly varying star 



formation histories with added nebular emission. To con- 
sider more carefully what this means, we need to understand 
whether these median age estimates describe the typical age 
of the galaxies, or whether they are driven by any other un- 
certainties. Also, we need to understand what the young 
ages from smoothly-varying SFHs reveal. 

8.2.1 Do any of the object SEDs require old stellar 
populations? 

To address the first consideration, we note that there is an 
inherent degeneracy in the two-component SFHs which pre- 
vent one from being able to distinguish between a fairly 
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Figure 10. The probability distribution functions (PDF) for the age of each object combining the results for the two-component SFHs 
with and without nebular emission (template sets D+E). 



Table 8. A table showing the f([OIII])/f(Hl3) ratio that would 
be required to produce the observed blue 3.6/^m — 4.5/im colours 
for four objects in our sample if the line emission were fueled by an 
AGN with X-ray luminosity at the sensitivity limit of the Chan- 
dra 4Ms data. The errors are determined from the photometric 
uncertainties in 3.6 fim and 4.5 fim fluxes. 



ID 


f ([Oil I}) /f {Hp) 


2 


28±13 


18 


11±6 


157 


55±21 


248 


57±28 



young (of order 100 Myr) burst that dominates the rest- 
frame optical flux, and an older burst that is mostly masked 
by the young star-forming component. The old median ages 
derived from the two component templates are therefore 
likely biased by this degeneracy and it is more meaningful 
the ask the question: Do we require old stellar populations 
to adequately describe the SEDs? 

It is clear from the object PDFs that most of the objects 



do not actually require stellar populations with particularly 
large ages. A few of the objects do show quite low probabil- 
ities of harbouring stellar populations much younger than 
100 Myr (e.g. 6, 8, 248), however. Objects 8 and 248, in par- 
ticular, consistently give fairly high minimum ages (lower 
age limit within la uncertainties) from fitting with the dif- 
ferent template sets. Plotting the object colours against the 
model colours (Figs [8] and [9| shows that the 3.6/im — 4.5/xm 
colour of object 248 is further than la from the template 
colours (the object with the bluest 3.6/^m — 4.5/xm colour). 
Given that the templates can try to reproduce red UV-to- 
optical colours by either boosting the nebular contribution 
or allowing for an old underlying population, it is possible 
that, if these objects have stronger [OIII] emission than ac- 
counted for in the models (as discussed in Section 8.1.1), 
the models would strive to reproduce this flux with an older 
underlying population. 
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8.2.2 Do any of the objects show evidence of being 
dominated by very young stellar populations? 

It is also interesting to note that three of the objects show 
distinctive peaks at extremely young ages (13, 26, 27) where 
the SEDs are well described by a burst with some recent star 
formation. With the current parameterisation, the peak for 
object 26 is likely to be unphysical given that the minimum 
Nl v c limit could not be applied (see Section 6). Object 13 
has a best-fitting model that requires reddening by dust and 
so the minimum Nl v c limit is likely to be underestimated. 
Each of these three PDFs also allow the full range of ages 
within la confidence. Raising the fraction of mass allowed 
within the star-forming component, however, would likely 
produce a high enough flux of Lyman-continuum photons to 
describe the Lya emission, and so it is important to check 
whether other objects might show such peaks at young ages 
using a constant SFR model with added nebular emission. 
We find slightly increased probability of very young ages (< 
50 Myrs) for many of the objects, although only objects 18, 
24 and 26 show a more significant peak at these young ages. 
However, these three objects have very uncertain derived 
parameters, either due to added degeneracies due to large 
amounts of reddening (24) or because the filters containing 
any Lya emission had to be excluded from the fit (18, 26). 

Using a template set consisting of only the smoothly- 
varying SFHs with added nebular emission, we find that two 
of the objects with very blue 3.6^tm - 4.5/^m colours prefer 
very young ages (~ 10 Myr) . This is reminiscent of the pos- 
sibly un-physically young ages (~ 3 Myr) fitted by Ono et al. 
(20101 to stacked observations of LAEs at z ~ 5.7. In fact, 



ages may no-longer be required. New evidence suggests that 
even Green Peas (low redshift systems showing evidence of 
a strong starburst and young ages) show older, underlying 
stellar components with low surface brightness that cannot 



four of the objects with the reddest J - H colours would also 
show young best-fitting ages but with broader uncertainties 
that extend to ~ 100 Myr, attributable to uncertainties in 
extinction. 

It is worth noting, however, that by insisting on a prior 
of smoothly varying SFHs, the only way to fit to the blue 
3.6/im — 4.5^im colours observed in a subset of these ob- 
jects would be with an extremely young population, as the 
nebular emission line equivalent widths fall quickly as the 
underlying population builds up. This would suggest an age 
dichotomy in the LBG population, with some of the ob- 
jects having just been formed with very low masses and 
very high sSFRs, and some at higher masses (those with- 
out blue 3.6/im — 4.5^im colours), still building up their 
mass. It is more likely, however, that the smoothly varying 
SFHs are primarily fitting to the dominant component by 
luminosity, the recent starburst. The young ages might be 
consistent with the expected duration of such a starburst, 
but the two component SFHs show that the SEDs also allow 
for an underlying population that can dominate the mass of 
the object. 

Although we cannot wholly discount the possibility that 
some of these objects have very young stellar populations 
with extreme [OIII] EWs (> 1000A), as observed in lower 
redshift galaxies such as Green Peas ( |Cardamone et aL 2009 1 
at 2 ~ 0.3 or the extreme emission line galaxies observed at 
z ~ 1.7 by van der Wei et al. (2011 1, we find that only two 
of the four objects with the bluest 3.6/xm - 4.5/xm colours 
show high probability of young ages when fit with smoothly- 
varying SFHs with added nebular emission. With a boosted 
/ ([OIII])/ f {Hp) ratio in the nebular models, such young 



be accounted for in SED fitting with a smooth SFH ( Amorin 
et al.|2012 1. 



8.2.3 Summary 

In summary, our estimates of the median best-fitting ages 
for our sample of z ~ 6 LBGs are of order 200-300 Myr, al- 
though the objects displaying the best age constraints give 
ages of ~ 50 — 200 Myr. We find no firm evidence for ex- 
tremely young stellar populations < 50 Myr and any objects 
showing old best-fitting ages (> 300 Myr), that also have 
colours that are well described by the template set, can also 
be described by stellar populations as young as ~ 100 Myr 
within the la confidence contours. 



We note that Schaerer & de Barros ( 2009 1 investigated 



the importance of nebular emission on the derived ages of the 
|Eyles et al.| ( |2007[ ) sample, of which four had spectroscopic 
redshifts at the time, and found that the best-fitting ages 
from templates with nebular emission gave a mean age of 
120 Myr (with some ages of order ~ 20 Myr) compared to 
500 Myr without nebular emission. Here, for a larger sample 
of spectroscopically confirmed objects with improved near- 
Infrared photometry from CANDELS, we do not see such a 
dramatic decrease in average ages between models without 
and models with nebular emission. This we attribute to the 
two component parameterisation of the SFH that allows for 
an underlying older stellar population that may dominate 
the mass, as discussed in the previous section. 



8.3 The Stellar mass-SFR relation 

Fig. 1 1 1 1 shows stellar mass vs. SFR for the full sample, de- 
rived using the full set of smoothly varying SFHs (a), as well 
as the const ant + burst templates without and with nebular 
emission ((b) and (c) respectively). In producing these plots, 
any models without current star formation (bursts) are dis- 
counted when determining the new best-fitting parameters 
and errors, which are the full range of masses and SFR of 
models with fits within Ax 2 = 1- 

The median sSFRs are plotted as the solid lines from 
the plot and are 1.7 Gyr" 1 , 1.8 Gyr -1 and 4.2 Gyr -1 for 
template sets C, D and E respectively. For template sets D 
and E this median excludes the objects without any Lya 
constraints, because the minimum SFRs are not restricted. 
The derived parameters from using template set D 
(b)) have much larger uncertainties in both mass 



11 



(Fig. 

and SFR, mainly introduced by the degeneracies between 
fitting the rest-frame UV light primarily by the burst com- 
ponent or the star-forming component. They do, however, 
illustrate that if the optical flux is dominated by low-mass, 
longer-lived, stars that also dominate the stellar mass, then 
the set of smoothly varying SFHs (i.e. Fig. |11^ ,) limit the 
lower possible values of SFR. This can contribute to im- 
posing a fairly well defined correlation when plotting sSFR 
against M*. If episodes of star formation occur over very 
short timescales, however, then the instantaneous SFR may 
not be closely coupled to the stellar mass on an object-by- 
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Figure 11. SFR vs. Mass derived using the different template sets as described in the captions. Errors shown are 68% confidence limits 
derived from A\ 2 = 1 from all models with SFR> M0yr — 1 . The three objects for which the SFR constraint from Lya was not applied 
(see Section 6) are plotted as triangles in panels (b) and (c). The line in each panel shows the median sSFR derived from that template 
set; 1.7 Gyr — 1 , 1.8 Gyr -1 and 4.2 Gyr -1 in (a), (b) and (c) respectively. The median sSFR for the two component models are derived 
excluding the three objects without Lya constraints (triangle symbols). 



object basis, even if the average properties of a population 
suggest otherwise. 

There is still the possibility, however, that red UV to 
optical colours are driven by emission-line contributions to 
the IRAC bands and not by the older stellar population. 
We see from Fig. 11 (c) that, compared to the results from 



the smoothly varying SFHs, there is greater scatter in the 
best-fitting parameters when nebular emission is included, 
as recently observed bylSchaerer & de Barros (20111. This 



is because there is a wider range of allowed masses, with the 
range of possible masses extending to lower values due to 
the uncertainty in the contribution of the nebular emission 
to the optical fluxes. 

We find that the estimated typical sSFR of the sample 
does increase using templates with added nebular emission, 
although by less than a factor of two. It is important to 
note, however, that our observations suggest that the ingre- 
dients of the nebular emission models may not be wholly 
representative of the sample, given that the colour space 
of the observations is not fully covered by the templates. 
A higher f([OI II])/ f(H f3) ratio would better represent the 
bluer IRAC colours in the sample, potentially leading to a 
higher mass and lower SFR estimate (because the colour 
relies on the EW of [OIII]). 



red line shows the relation derived in Gonzalez et al. (2011 1, 
also from an LBG sample at z ~ 4. 



The errors in mass are taken from the range of ac- 
ceptable mass within A\ 2 = T Our results agree qual- 
itatively with the flattening of the slope evident in the 
|Stark et al. | ( |2009[ ) relation, while at magnitudes fainter 
than M1500 = —21.5, the relationships from both papers 
agree. Any flattening to high UV luminosities could not be 
observed in the |Gonzalez et al.| ( |20lij ) results, however, be- 
cause their template set was restricted to a constant SFR 
model. 



The masses derived from the templates that include 
nebular emission extend to lower values than estimates from 
the other two template sets, although they still quantita- 
tively agree with the z ~ 4 relationship. We do, however, 
find a marked increase in derived sSFRs from fitting with 
template set E compared to the results derived using tem- 
plate set C (Section 7.2) that is not so readily observed in 
the M* — M1500 relation, especially when full uncertainties 
in mass are considered. 



8.4 The stellar mass-UV luminosity relation 

We have seen that the errors involved in determining the 
SFR for these galaxies are large given the degeneracies be- 



From the apparent lack of evolutionin the M» — M15 



tween SFR and extinction. Stark et al. (2009) chose to in- 



vestigate any evidence of mass growth for successive pop- 
ulations of LBGs between 4 < z < 6 via the M* — M1500 
relation in order to distinguish between various scenarios 
of LBG evolution. This can also reveal whether the typical 
sSFR and stellar masses change significantly over this time. 

Fig. [12] shows the best-fitting masses derived from tem- 
plate sets C, D and E (panels (a), (b) and (c) respectively) 
vs. UV absolute magnitudes for this sample. Over-plotted in 
green is the median stellar mass-UV relation derived from 
the Stark et al. (20091 z ~ 4 sample (from masses derived 
using an exponentially decreasing SFH with r = 100 Myr 
and corrected to the Chabrier IMF; see their Fig. 9) and the 



relation, Stark et al. ( 2009 1 argue that it rules out a steady 



growth scenario in which the LBGs from the higher redshift 
populations continue to grow at a constant SFR, as the de- 
rived ages imply that they would populate successive gener- 
ations, shifting the normalisation of the lower redshift pop- 
ulations to higher masses at constant UV luminosity. They 
suggest instead that short duty cycles implied by exponen- 
tially decreasing SFHs fit the observations better, although 
this would mean that each of the objects observed would 
have had a higher SFR in the past. Our results also display 
little evolution between our z ~ 6 sample and the Stark 
et al. (20091 z ~ 4 sample, consistent with the theory that 



LBGs selected at each redshift are an independent popula- 
tion. We also see, however, that simple two component SFHs 
describe the SEDs well and note that more stochastic SFHs 
would also fulfill the requirement for short duty-cycles. 
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Figure 12. Mass vs. UV absolute magnitude derived using different template sets (as indicated by the panel captions). Errors shown 
are 68% confidence limits derived from A% 2 = f. Over-plotted in green is the relation derived from the z ~ 4 sample of Stark et al. 
(2009) and the red line shows the relation derived at z ~ 4 in Gonzalez et al. (2011) 



9 CONCLUSIONS 

By performing SED fitting to a sample of spectroscopically 
confirmed LBGs, we are able to remove the uncertainty in 
the redshift from the derived parameters, as well as the 
contribution of Lya to the observed photometry in all but 
two cases. Moreover, deep near-IR imaging supplied by the 
CANDELS survey and deep IRAC imaging from Spitzer 
SEDS has allowed us to accurately constrain the rest-frame 
UV to optical photometry. Fitting with templates cover- 
ing a wide range of SFHs, metallicity, extinction, as well 
as whether or not they include nebular emission, we find 
that: 

(i) When fitting with a constricted template set with a 
single metallicity and no extinction, the derived best-fitting 
ages show 3/13 of the objects with ages > 500 Myr. Expand- 
ing the template set to contain a wide range of SFHs, metal- 
licities and extinction, however, shows large uncertainties in 
the derived ages, although the derived masses are shown to 
be robust to within a factor of two, in accordance with other 
studies (e.g. |Yabe et al.|2009||Pacifici et al.|2012 l. 

(ii) Motivated by an observed tension between the ob- 
served photometry and models that caused objects with 
strong apparent Balmer-breaks but blue UV slopes to be fit 
with extremely old templates, simple two-component mod- 
els were used to decouple the assembled stellar mass from 
the instantaneous SFR. Fitting these templates, both with 
and without nebular continuum+line emission, shows that 
6 of the 13 objects are marginally better fit without nebu- 
lar emission, while the remaining objects show better fits to 
their blue 3.6pim - 4.5 /jm colours when nebular emission is 
included. The bluest 3.6/im - 4.5 /im colours are not well de- 
scribed by any of the models, although a moderate increase 
in the flux of [OIII] compared to H/3, consistent with mea- 
sured flux ratios in some z ~ 3 LBGs, would better describe 
the observed colours. Masses derived using templates with 
nebular emission are systematically lower than those from 
templates without nebular emission, with a median offset of 
0.18 dex. 

(iii) Full analysis of the individual probability distribu- 
tion functions demonstrates the poor constraints that can 
be placed on stellar population age on an object by object 
basis. Without spectroscopic redshifts and the extra infor- 
mation provided by the Lya flux these constraints would 



be poorer still. The combined probability distributions for 
the two-component models, both with and without nebular 
emission, show no firm evidence for extremely young stellar 
populations (< 50 Myr) and no objects that require popu- 
lations with ages > 300 Myr to describe their SEDs. The 
three objects with the strongest constraints give ages in the 
range ~ 50 - 200 Myr. 

(iv) Results from all of the different template sets show 
that specific SFRs derived using templates with nebular 
emission are moderately higher than those derived without 
nebular emission (median sSFRs of 4.2 Gyr -1 compared to 
1.8 Gyr" 1 respectively), although the uncertainties are large 
and the scatter is greater for the models including nebular 
emission, in agreement with recent results by [Schacrcr & 
de Barrosl (120111) andlde Barros et al.l (12011 1. 



(v) The derived M* — Luv relation displays little evidence 
of evolution in the LBG population between 4 < z < 6 
which can naturally be described either by a more stochastic 
parameterisation of the SFH, with short durations of star- 
formation, or an exponentially increasing SFH. 
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APPENDIX A: PHOTOMETRY 

In this appendix we provide the multi-wavelength photome- 
try for each of the objects in our sample. In Table Al we list 
the photometry for the eleven objects within the GOODS-S 
field and the photometry for the two objects within the UDS 
is listed in Table A2. 
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Table Al. Photometry of the objects in GOODS-S field. Column 1 gives the object ID while columns 2-10 give the AB magnitudes and 
corresponding errors in the following filters: the HST ACS B435, V6O61 *775 ar >d Z850 _Dan d filters, HST WFC3 J125 and Hieo filters, the 
HAWK-I Kg— band and IRAC 3.6^tm and 4.5^tm band filters. For the purposes of the SED fitting the photometric errors for the optical 
and near-infrared bands have been allocated a minimum error of 10%, and the IRAC fluxes have been allocated a minimum error of 
20%. For undetected objects in any of the bands, the 2a limiting magnitudes are reported. All magnitudes have been aperture corrected 
to the same enclosed flux provided by a 0.6"-diameter circular aperture at the spatial resolution of the HST WFC3 Higo imaging. The 
aperture corrections needed to convert the photometry into total magnitudes (based on the -H16O measurements) are listed in column 11. 



ID 


S435 


V606 


*775 


2850 


J125 


Hi 60 


K s 


IRAC1 


IRAC2 


aperture 
correction 


2 


>28.44 


>28.68 


27.38±0.30 


26.04±0.10 


25.89±0.10 


26.02±0.10 


n/ a 


25.20±0.20 


25.66±0.25 


0.60 


5 


>28.43 


>28.37 


26.41±0.12 


24.78±0.10 


24.79±0.10 


24.90±0.10 


n/ a 


24.38±0.20 


24.57±0.20 


0.20 


6 


>28.53 


>28.78 


>28.15 


26.18±0.10 


26.50±0.10 


26.62±0.13 


27.22±0.94 


25.90±0.20 


25.37±0.20 


0.32 


8 


>28.65 


>29.13 


26.83±0.15 


25.43±0.10 


25.45±0.10 


25.53±0.10 


25.49±0.29 


24.63±0.20 


24.77±0.20 


0.31 


13 


>28.62 


>28.80 


>27.53 


26.63±0.16 


26.53±0.13 


26.14±0.10 


26.22±0.30 


25.72±0.20 


25.68±0.20 


0.25 


18 


>28.84 


>28.65 


>28.29 


26.30±0.10 


26.94±0.18 


26.96±0.26 


27.48±1.18 


26.36±0.20 


26.90±0.47 


0.20 


23 


>28.08 


>28.51 


>27.83 


26.37±0.14 


25.54±0.10 


25.69±0.10 


n / a. 


24.82±0.20 


24.75±0.20 


0.48 


24 


>28.54 


>28.68 


27.01±0.19 


25.87±0.10 


25.37±0.10 


25.26±0.10 


24.92±0.22 


24.17±0.20 


23.68±0.20 


0.49 


25 


>28.52 


>28.61 


>27.82 


26.23±0.14 


25.98±0.10 


26.11±0.10 


26.15±0.44 


25.30±0.20 


25.32±0.20 


0.73 


26 


>28.63 


>28.76 


>27.96 


26.55±0.15 


26.48±0.10 


26.29±0.10 


n / a. 


25.70±0.37 


25.82±0.41 


0.24 


27 


>28.10 


>29.16 


>27.89 


26.14±0.10 


26.22±0.10 


26.22±0.10 


26.22±0.54 


25.69±0.20 


25.86±0.20 


0.41 



Table A2. Photometry of the objects in UDS field. Column 1 gives the object ID while columns 2-11 give the AB magnitudes and 
corresponding errors in the following filters: the Subaru B, V, R, i and z-band filters, HST WFC3 J125 and Hieo filters, the UDS WFCAM 
if— band and IRAC 3.6/^m and 4.5/^m band filters. For the purposes of the SED fitting the photometric errors for the optical and near- 
infrared bands have been allocated a minimum error of 10%, and the IRAC fluxes have been allocated a minimum error of 20%. For 
undetected objects in any of the bands, the 2tr limiting magnitudes are reported. All magnitudes have been aperture corrected to the 
same enclosed flux provided by a 0.6"-diameter circular aperture at the spatial resolution of the HST WFC3 H160 imaging. The aperture 
corrections needed to convert the photometry into total magnitudes (based on the Hiao measurements) are listed in column 12. 



ID 


B 


V 


R 


i 


z 


J125 


-H16O 


K 


IRAC1 


IRAC2 


aperture 
correction 


157 

248 


>28.78 
>28.78 


>28.00 
>28.00 


>28.31 
>28.31 


>28.24 
>28.04 


25.05±0.10 
25.33±0.10 


25.31±0.10 
25.47±0.10 


25.55±0.10 
25.79±0.10 


25.63±0.53 
25.85±0.60 


24.69±0.20 
24.46±0.20 


25.52±0.28 
24.94±0.27 


0.21 

0.20 
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